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Abstract: Disorder has long been a difficult subject in condensed matter systems and tfie 
The replica method is a well-known tool in this field. Implementing the replica method the 
AdS/CFT correspondence has been proposed and discussed in literatures. We point out, for 
any CFT that has a holographic dual and to the leading order of the large-A*" expansion, 
the corrections due to the presence of random disorder to any connected correlation functions 
vanish identically, provided that the disorder strength is normalized as discussed in literatures 
and that the symmetry among replicas is unbroken. Same must hold true to any observables 
that are determined by the connected correlation functions through a linear relation. This 
behavior resembles strongly that of a free theory where disorder is coupled to the fundamental 
field. We demonstrate this by both the means of holographic principle and field theory 
analysis in a toy model. We also propose ways of evaluating the non-zero sub-leading effects 
perturbatively in terms of the disorder strength and discuss a novel possibility of defining a 
new holographic dual if we adopt a different normalization for the disorder strength. 
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1 Introduction 

Random disorder is common in many-body systems. In condensed matters, spatial inhomo- 
geneity, or impurities, is almost always present. In certain circumstances, they are the key 
element responsible for important physical properties. A famous example is the DC conduc- 
tivity, which is infinite theoretically in a translational symmetric system and only becomes 
finite because of the presence of impurities in real material. Despite its importance, disor- 
der in quantum system remains a difficult problem, particularly when the system is strongly 
coupled. 

One often treats the impurities statistically. It is first introduced into the theory as a 
classical source coupled to a certain operator in a QFT, and then the system is averaged 
over a probability distribution functional for the random source. There are two main diffi- 
culties in this process. First of all, systems deformed by an arbitrary source loose the spatial 
translational symmetry and become difficult to handle, and secondly, taking the average is 
not trivial. Two well-known methods to overcome both difficulties are the so-called "replica 
method" and the "Grassmann field method" [1]. Both aim at translating the said problem 
into an ordinary QFT calculation such that more tools are at one's hand. The replica method 
involves duplicating the theory into multiple copies, introducing a mixing term that couples 
them, and taking the limit that the number of copies approaches zero. While the replica 
method is an elegant formalism and a vast amount of literatures are devoted to the topic, it 
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is not always easily implemented in many physically interesting systems, partially due to the 
bizarre limit one has to take. 

Ever since its discovery [2-4], AdS/CFT correspondence has been applied widely to many 
types of systems. It did not take long before people attempted to utilize its power to advance 
our understanding in condensed matter theories (see, for example, [5-8] and the references 
therein). Given that fairly simple prescription for holographic duals for CFTs with multi- 
trace deformations was known [9] , for reasons to become clear below, it naturally leads to the 
proposal of implementing the replica method in holography and the hope for making progress 
in disorder problems [10]. Some preliminary attempts based on this idea to evaluate the DC 
conductivity have been made in [11]. 

Apart from all the technical subtleties, part of the weakness of the discussions in [11] is 
that the bulk geometry considered there was pure AdS space, corresponding to a system at 
zero temperature and with a zero charge density. One would much prefer to consider models 
that are more realistic, in which both the temperature and the charge density are finite. 
Indeed, holographic models of such kind exist. One of the famous stories in this category is 
the model for holographic superconductors where a hairy black hole replaces the pure AdS 
geometry [12-14]. Chances appear to be much better that one finds more interesting results 
there since it is possible to have a non-zero background profile for some bulk scalar field, 
and if one lets the random disorder couple to the operator dual to that field, the background 
geometry as well as all the perturbations above that are sensitive to the presence of disorder, 
leading to direct corrections to boundary-to-boundary correlation functions. 

In this note, we show that, despite the above naive expectation, random disorder has zero 
corrections toward any connected correlation functions in the leading order of the large-N 
expansion, provided that the symmetry among replicas is not broken. The last condition is 
a subtle assumption and systems known to violate it exist [15]. Those must be studied in 
separation. It must be emphasized that observables depending on the disorder in a non-trivial 
way certainly exist (such as that in [16]) and may be evaluated using the replica method. But 
the connected correlation functions themselves in the leading order of N is independent of 
the disorders that have a Gaussian distribution profile, and any quantities determined by the 
connected correlation functions through a linear relation, such as the DC conductivity, are 
not corrected either. These results resemble very strongly the case where a single field is 
coupled to a random external source in a free QFT, even though CFTs with semi-classical 
AdS duals are strongly coupled. We explain our statement using the replica method both by 
the holographic duality as suggested in [10] and field theory analysis perturbatively in terms 
of the strength of random disorder for a toy model that has a similar large-A^ expansion. 
Both methods agree well. 

At the sub-leading order, random disorder may lead to non-trivial corrections. Based 
on our understanding gained from the field theory model, we propose ways to evaluate those 
effects holographically. Unlike the leading order effects, however, we are only able to do so 
perturbatively in terms of the strength of disorder. When the disorder strength is normal- 
ized in the standard way as discussed in literatures, higher order contributions are further 
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suppressed by more powers of and become less and less important. In fact, it is not use- 
ful to evaluate them if one does not also include many other competing effects, such as the 
non-planar contributions at the same time. 

We also discuss the possibility of making the disorder stronger as N ^ oo such that 
the non-vanishing corrections can compete against other leading order physics in the large- 
expansion. In fact, it appears to be the most natural and interesting choice since higher order 
contributions in terms of the disorder strength are no longer suppressed and all appear at the 
leading order in the large- expansion. We will show that this choice amounts to an unusual 
normalization for the double-trace deformation in a CFT. Such a normalization is usually 
considered pathological because the deformation destroys the large- expansion in ordinary 
CFTs. But in the context of the replica method, the sickness is cured by the appearance of a 
new small parameter, namely the number of the replicas. In the queer limit that this number 
approaches 0, a consistent large-N expansion may exist, leading to the hope that a new 
holographic model can be defined and it captures the effects of disorder nonperturbatively. 
Unfortunately, such a theory remains unknown to us. 

The rest of the paper is organized as the following: in the next section, we briefly ex- 
plain the replica method emphasizing the "magic limit" in which it "converts" disconnected 
correlation functions into properly normalized connected ones. In Sec. 2.2 we implement the 
replica method in holography as it is proposed in [10]. In Sec. 2.3, we present the leading 
corrections caused by the disorder to the one- and two-point correlation functions using the 
AdS / CFT correspondence and show that they vanish when the symmetry among the replicas 
are preserved. In Sec. 3, we reproduce the same results perturbatively in the matrix model 
by field theory methods. In Sec. 4, we discuss how one may evaluate the sub-leading non-zero 
effects, and the possibility of making them stronger and defining a new type of holographic 
dual for CFTs with an unconventionally normalized double-trace deformation. We make some 
further discussions in Sec. 5 and present a few more details in the appendix. 

2 The old holographic replica method and the correlation functions 

In this section, we briefly introduce the replica method and provide an essential set of identities 
needed for later discussions. We also explain how one may implement the replica method 
in AdS/CFT as it is discussed in literatures, where the disorder strength is normalized most 
conveniently such that its leading contributions can be easily evaluated using known methods. 
Unfortunately, as we show in the last part of this section, if one is only interested in connected 
correlation functions, such contributions, even though being non-perturbative in terms of the 
disorder strength, are simply zero, assuming that the symmetries among replicas is preserved. 

2.1 A brief introduction to disorder problem and the replica method 

The physics problem to consider is evaluating the correlation functions in some QFT at 
the presence of a random external source averaged over its probability distribution. Let us 
denote the random source as V{x) and couple it to an operator 0{x) in the QFT. For any 
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fixed V{x), the action for the QFT is deformed into Sy = S[X] + J d'^xV{x)0{x), where S[X] 
is the action for the original theory and X represents collectively all the fundamental fields. 
As an example, the averaged expectation value for O is given by 



(0(x)) = J B[V]P[V 



Here, P[V] is a predetermined probability distribution functional for V{x). In real materials, 
such a distribution functional is usually localized at the time-independent configurations for 
V{x), corresponding to impurities that only break spatial homogeneity but preserve the energy 
conservation law. Throughout this discussion, we would choose to ignore this subtlety and 
simply assume that P[V] takes the simplest possible Gaussian form as 
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where we introduced the "strength of the disorder" /, assumed to be a constant below. To 
make contact with more realistic physical systems, one should promote / to an operator, or 
in the momentum space, a general function of the frequency and momentum. In particular, 
one may multiply / by 6{uj), where uj = is the frequency, so that V{x) is restricted to 
be time-independent. In the real-time formalism, this delta-function is needed to recover 
unitarity. 

As usual, we define the generating functional 

^S[X]+f d'^x{0{x)V{x)+0{x)J(x)) 



and write 



Zv[J] = Jb[X] 

J 6J(^x) 



{0{x) 

The replica method is invented based on the following identity 

n->0 n 5J[xi)5J{x2) ■ ■ ■ 5J{Xm) 5J{xi)5J{x2) ■ ■ ■ 5J{Xm) 
Z\jy^ can be expressed explicitly by rewriting the same path-integral n times: 

Here, we duplicated the same theory, including both the fields and the action, for n times, 
and hence the name of "replica method" . While n is defined as an integer in this procedure, 
we must assume that the physical observables can be understood as analytic functions of n 
so that the limit n — t- is defined. 

Granted that the limit in Eq. (2.1) makes sense, the magic of replica method 

^ I n n n \ ° 

hm -( VO,(xi) V0i(:E2)---V0.(x„) ) = ( 0(xi)0(x2) . . . O(x^) ^ (2.2) 

\ i=\ 2=1 j=l / 
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follows immediately. The left-hand side of this equation contains the unnormalized discon- 
nected correlation functions, denoted by the notation ( * )° throughout this manuscript, eval- 
uated in the theory with replicas, and the right-hand side is the normalized totally connected 
correlation functions, denoted by the notation ( * )'^, in the original single-copied theory. The 
former, if divided by n, approaches the latter in the limit n — )• 0. 
Averaging both sides over the distribution we find 

(0(xi)0(x.)...0(x„)r=limi / D[y]e-^^'^''-^(-)^l ^^^ J7\^'^K u 

n^onj n6J{xi)6J{x2) . . .6J{Xrn) 

Assuming that one can interchange the order of the limit n — )• and the functional integral 
for V, we define 



where 



is the quadratic deformation, or in a CFT, the "double-trace" deformation, and arrive at our 
master formula 

1 



{0{xi)0{x2)...0{Xrr,)y = lim 



n-s>0 n 5J{xi)5J{x2) ■ ■ ■ 6J{Xm) 

\ i j k I J 

We add a subscript "/" for the last correlation function in the above equation to emphasize 
that it is evaluated in the theory with replicas plus the quadratic deformation (2.3). 
Particularly, for the vacuum expectation values, we have ^ 



(0) = lim-/ Vo\ = lim(Oi)° , 



/ 

and the connected two-point correlation functions 



I Th TL \ 

{0{x{)0{x2)r=Y^^J-{Y,0,{x,)Y,0,{x2) ) . (2.4) 



i=l j=l 

It is easily verified that 

( 0{xi)0{x2) ) = lim ( Oi(xi)Oi(x2) )) 
( 0{xi) ) { 0{X2) ) = lim ( Oi{xi)02{x2) )) 



(2.5) 



'^The last step assumed the symmetry among replicas are unbroken. 
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The particular values of the two subscripts for both operator Oi on the right-hand side of 
both equations above are not important other than that they must equal for the first equation 
and different for the second. 

It is sometimes useful to make the follow rotation among the replicas: 



Oi 



(2.6) 



where a 



nj 



-4= and 



1, l<j<i 
-i, j = i + l 

0, j > i + 1 



for i < n , 



easily verified to be orthogonal. Straightforward calculations show 

1 r 



{OnOn)° 
{ OnOn-1 )° 



n 

1 r 

n 



OnOn) + {n - 1) 
OnOn ) — { On-lOn-1 



By Eq. (2.5), in the limit n — )• 0, we have 



{00) = lim ( OnOn )° = lim ( 0„_iO„_i ) + lim - 

n^O n^O \ I n^O n 



[OnOn 



y On-lOn-1 



{O) {O) = Yim {OnOn-l)° 

71— !>0 



lim — 

n-s>0 n 



OnOn ) - ( On-lOn-l 



(2.7) 



Therefore, we also find ( OO )^ = limn_j.o On-iOn-i j ■ Comparing this to Eq. (2.4) implies 
the replica method is only self-consistent if the following two limits coincide: 



lim ( OnOn 
n-s>0 



lim ( On^lOn-l 
n->0 



In fact, it is required that the difference between ( OnOn ) and ( On-iOn-i ) is 0{n) so 



that ( O ) ( O ) would not diverge, as shown in Eq. (2.7). 

We mention that rotating into the basis of Oi is not extremely useful in general since both 
the action S and the operator 0[X] are typically nonlinear functionals of X, and not only it 
is non-trivial to implement such a rotation in terms of the fundamental fields Xi but also it 
creates more mixing in other parts of the action leading to no simplifications. Only in the 
trivial case explained in the Appendix A, it is helpful simply because the action is quadratic. 
It turns out, to the leading order of the large- expansion, a class of problems that can be 
treated by holographic methods discussed below appears to present another example where 
such a rotation is useful, leading to just as trivial results. 
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2.2 The replica method in AdS/CFT 

While the replica method is a beautiful idea, it is difficult to implement in a generic QFT, 
particularly when it is strongly coupled. On the other hand, to the leading order in the 
large- expansion, holographic dual for a CFT deformed by a multi-trace term was known 
[9] and a vast amount of literatures were devoted to studying those theories [17-25]. One 
naturally hopes that implementing the replica method in the AdS/CFT correspondence may 
be tractable and bring us to better understandings to disorder problems. 

First, let us set the number of the replicas to 1 and briefly mention how the multi-trace 
deformation of the form 

bL = W[0] (2.8) 

is treated holographically. Here O is some single-trace operator in the CFT and W[x] is a 
polynomial of x. 

The bulk geometry we have in mind is an d -|- 1 dimensional asymptotically AdS space. 
It may be an AdS charged black hole, or an hairy charged black hole as discussed in [12-14]. 
We will not need to specify the geometry in details in this discussion. 

Let ip be the bulk field dual to O. in such an asymptotic AdS space. Near the boundary, 
if the metric approaches the form ds^ = z~'^{dz'^ + dx'^dx^), ipi behave as 

where A± = d/2 i: y^d'^/A + m?. We identify A_ with the conformal dimension of O. By 
the holographic principle, when -0 satisfies the boundary condition a = near z = 0, /? is 
identified with ( O ) in a CFT without deformations. When the CFT is deformed by a multi- 
trace term (2.8), it was proposed that, to the leading order of A, tp must satisfy the modified 
boundary condition [9] 

5W[I3] 

and the rest of the AdS/CFT dictionary remains unchanged. The double-trace deformation 
usually breaks the conformality and lead to non-trivial RG flow [9, 25]. We have above 
restricted ourselves in the window such that the deformation is relevant, i.e. 2Ao = 2A_ < d. 
Correspondingly, the mass for the bulk scalar must be within the window — ci^/4 < rm? < 
— (i^/4 -|- 1 so that alternative AdS quantizations are allowed [26] 

Given this simple prescription, it is suggested that the replica method may be imple- 
mented in the AdS/CFT correspondence as the following [10]. Duplicate all the fields in the 
bulk, including the metrics, by n times and set the boundary condition for the fields dual to 
the operator O as mentioned above. For example, consider a phenomenological model given 
by the Lagrangian 

n ^ 

C = y] [R{gi) + 6 - -F^.F^'^Ai) - {d^^i - iqA^i^iXd''^ + iqAf^^) 

i ^ (2.9) 

-mViV'I -A|Vi|^ + ...] , 
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where we have set the bulk curvature and the Planck mass to 1, and introduced n copies for 
every field, including the metrics gi, the gauge fields Ai = A^dx'^, and the complex scalars tpi. 
The scalar would be dual to the operator O in the CFT that is coupled to the impurity ^. 
Toward the boundary, ^/^j ~ aiZ^+ + f3iZ^- , for which we should set the boundary conditions 
as 

n 

ai = fY,^i- (2-10) 

i 

Such boundary conditions not only mix different copies of the scalar fields, but, through loop 
corrections, also mix all the gauge and gravitational fields, making all but one linear combi- 
nation of the replicas massive, reflecting the fact that only one conserved energy-momentum 
tensor and U{\) current are preserved on the CFT side due to the mixing double-trace defor- 
mation. We will evaluate the correlation functions in this theory by the standard AdS/CFT 
duality and take the limit n — t- 0. 

We should mention discussions on applying holographic models for CFTs with double- 
trace deformations to condensed matter problems exist in literatures, and they are not always 
motivated by the replica method explained above (see, for example, [27]). When the multi- 
trace deformations are introduced for other reasons, one is probably not interested in taking 
the limit n — )• 0, or such a limit simply does not exist, in which case the following discussions 
of course do not apply. 

2.3 The vanishing leading order corrections 

Let us first examine the possible corrections to the vacuum expectation values. For this pur- 
pose, we look for the background solutions to the equations of motion led by the Lagrangian 
(2.9). Schematically, we can write the equations as 

*,[V'f,5f,^f] = 0, E,[V'f,5f,^f]=0, M,[V'f,5f,Af]=0, i = l,2,...,n, 

where 'J', E, and M represent the Klein-Gordon equations, the Einstein's equations, and the 
Maxwell's equations respectively. The solutions for the metrics must be asymptotically AdS 
but may develop a black hole in the bulk. These equations are nonlinear and known analytic 
solutions are rare, but can often be solved numerically. These equations do not mix fields for 
different replicas, but the boundary conditions (2.10) do. 

We assume that the background solutions respect the translational symmetry in the d 
flat dimensions and are x^-independent. Given such assumption, all the equations above 
form a system of ordinary differential equations and fields in different replica share identical 
equations as well as boundary conditions. If the vacuum is unique, we must find the solutions 
for all replicas are identical. In the scenario that there exist multiple vacua in the bulk, the 
solution that corresponds to the least action should be favored and the action in the bulk 

■^Here, the operator O is charged, and so must be V that's coupled to it. The double trace deformation term 
would then become ^ (^^ O;)^ (Sj Oj^. Similarly, the correlation functions to evaluate should be (J O^O Y • 
For the simplicity of notations, we would ignore the complex conjugation though. 
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is simply a sum of those of all replicas, it is natural to assume that each replica picks the 
same solution. However, we should emphasize that there exist possible scenarios where this 
symmetry is broken if more complicated potentials are included, and when such phenomenon 
happens, the problems is substantially more complicated [10, 15]. We will not further discuss 
this possibility throughout this paper. Given these assumptions, mixing boundary conditions 
(2.10) automatically decouple and become equivalent to 

ai = nfPi, i = l,2,...,n. (2.11) 

It is manifest from Eq. (2.11) that the background profiles depend on the strength of the 
double trace deformation / only through the combination of / = nf. Hence, taking n — >• 
limit becomes equivalent to turning off the disorder. For example, ( O ) is given by 

(Oy = lim - / V Oi ) = lim ( Oi ) . = lim . 

\i=l If 

The limit n — t- and / — >■ is indistinguishable in the last expression and it is necessarily 
true that 

XW)) f = {o[x)) , 

as if there is no random disorder at all. Here, a key requirement is that the background 
solutions {ip'^\ g'^'' , A'^^) vary with respect to / smoothly as / approaches 0. It certainly holds 
true in the present case with the assumption that the solutions for all replica are identical, 
and in the limit / — )• 0, the standard boundary condition is recovered from Eq. (2.11). 

Now, let us turn to the two-point correlation functions. To this end, we should consider 
the fluctuations above the background profiles: gi — t- gf + 5gi, Ai — )• Af + 6Ai, and ipi — )• 
ipf + dipi, and solve for them by perturbing the equations of motions to the first order. 
Schematically, we can write: 

K[g'\A'\r']SXi = 0, (2.12) 

where 5X denotes the perturbations {5g, 6 A, 5tp) collectively, which may be thought of a 
column vector and X is a linear second order differential operator obtained by expanding the 
Einstein's equation, the Maxwell equation and the Klein-Gordon equation to the first order 
of the perturbations above the background {g'^^ , A'^^ ^ip'^^). Fields in different replicas must all 
be allowed to fluctuate freely and their perturbations are generically different. However, here 
we have the extra advantage that the perturbed equations Eq. (2.12) form a linear system. 
We may apply the rotation introduced in Sec. 2.1, quite similar to the way we solve the free 
theories outlined in Appendix A. Define, for each Xi S {^j, ^j, ■i/'i}) 

n 

Xi = ^^ aijXj , 
i=i 

where Ojj are given in Sec. 2.1. According to Eq. (2.4) and the AdS/CFT correspondence, 
we find 

{Oxix)Oxiy))cFT = lim ( Ox„ix)dxAy) )° = Ihn H^^^Jx, y) , (2.13) 
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where H^^j^^ is the boundary-to-boundary propagator for the field X„ and Ox denotes the 
operator that is dual to X. Some objections may be raised against the last step above. 
In the intermediate step, one is supposed to evaluate the unnormalized disconnected two 
point function for the operator Ox„) but by the AdS/CFT dictionary, the boundary-to- 
boundary propagator is identified with the normalized connected correlation function. In this 
context, the difference is not important. First of all, the normalization is irrelevant since 
a key assumption for replica method to work is that Z„ — t- 1 in the limit n — t- 0, where 
Zn is the partition function. While this limit is subtle as Zn is usually divergent and one 
must properly regularize it, provided the proper regularization is found, we may interchange 
the normalized and unnormalized correlation functions in the limit n — t- 0. Secondly, the 
disconnected correlation function differs from the connected one by a possibly nonzero term 

^ ^ • But this term is always proportional to n, as is easily verified, and vanishes in the 
limit n — )■ 0. A quick way to see this must be true is to recall that 

hm ( Ox„Ox„ ) = lim ( Ox^.^Ox^., ) 

is essential in replica method, and by construction (| Ox„_i ^ = Therefore, we may as 

well replace the field by Xn-i, which leads us to the conclusion below more quickly. 

In any case, the boundary-to-boundary propagator for Xn is given by the solution for 
6Xn- In the tilde-basis, the boundary condition is "unmixed". Depending on the correlation 
functions oiiG is intcrGstcd. in, bouiidcLry conditions for fields other tlian ifj^i are standard ones 
and only that for S'lpn is modified. For example, if one wish to evaluate Ox„ix)Oxn{y) ^ ; 
corresponding to the response of the system with respect to an infinitesimal change of the 
source that's coupled to O, one may set 

5an = 1 + nf6f3n , 

for 6^n, and make sure other fiuctuations vanish near the AdS boundary, after proper regu- 
larizations if necessary. The averaged two-point correlation function for the order parameter, 
{0{x)0{y) Y, is then given by lim„_j.o 5/3n- Both the boundary condition and the operator 
K depend on /, and K acquires its dependence through the background solutions. However, 
the appearance of / here in the boundary condition is again in the combination of / = n/, 
and taking the limit n — )• is still equivalent to taking the limit / — )• 0. The same is true 
for K simply because it happens to the background solutions as discussed. Hence, taking the 
limit n — )• completely eliminates the /-dependence in Eq. (2.13). 

In fact, as is alluded to already, the consistency of the replica method demands that 
the same limit is achieved by considering 5/3„_i instead, for which the standard boundary 
condition applies. The only /-dependence of 6(3n-i comes through K. We, therefore, are led 
to the conclusion even faster that the connected two-point correlation function are insensitive 
to the presence of the disorder, when the double-trace deformations is taken into account by 

^Once again, the symmetry among all replicas is used 
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the generalized boundary conditions as explained and with the assumption that the symmetry 
among the replicas is unbroken. Despite that it may appear counter-intuitive, this conclusion 
is explicitly verified by the results presented in [10]. 

We should mention that it is certainly not true that no observables depend on / exist. 
For example, both 

5f3n + in-l)5Pn-i] , 

6i3n - 5i3n-l 

depend on / non-trivially because the difference between (5/3„ and S/3n-i is 0{n). 

Similarly, to evaluate the current-current correlation function, we should solve for SAn 
with the boundary conditions 

6An = 1 , 6an = nfdPn , 

and all other perturbations approach near the AdS boundary. By the same reasoning, in 
the limit n — )• 0, the averaged two-point function should just become trivial: 

( J^(x) J^(y) = { Jf,{x)Ju{y) )}=o ' 

as if / = in the first place. 

So, given the said assumptions, we find disorder has no effect on the averaged connected 
correlation functions in the leading order of the large- expansion in any CFT that has a 
holographic dual described by (2.9). The way we arrive at this conclusion is very similar to 
how we solve the same problem in a free theory deformed by disorder as shown in Appendix 
A. Certainly, the symmetry among the replicas play an essential role which allows one to 
diagonalize the system in its holographic description very easily. 

3 A toy field theory model 

What we found above is not immediately intuitive considering the CFT is supposedly strongly 
coupled. We seek a field theory understanding toward the same problem in this section. 
Besides the very much needed reassurance to verify that the holographic principle does not 
break down, a field theory understanding might also allow us to probe further the sub-leading 
corrections, if any, not captured by the treatment presented above. 

It is not easy to reproduce the results given above in an arbitrary CFT deformed by 
the disorder. But, perturbatively in terms of the disorder strength /, we may make some 
attempts in a toy model that share some important properties with theories that have a 
holographic dual. We emphasize that when we claim to analyze the system perturbatively, 
we always mean that we do so perturbatively with respect to the disorder strength /, and 
make no assumption about the theory itself being strongly or weakly coupled or whether the 
theory can be treated perturbatively with respect to any other couplings. As long as we can 



(0(x)0(y))cpT = lim- 



( 0{x) )cFT 0(y)cFT = - 
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tune the parameter / to be sufficiently small, we assume that a perturbative expansion for 
the double-trace deformation, or more generally the multi-trace deformation, in the orders of 
/ can be defined even if the theory is strongly coupled. 

The toy model we present here is the matrix theory. While it is not a CFT, we may 
proceed since the discussions below solely rely on the assumptions that the disorder can be 
treated perturbatively and the theory has a large- A'^ expansion and well-defined t' Hooft limit 
in which correlation functions have the similar factorization properties as those in the matrix 
model. 

The theory is described by a Lagrangian: 

£ = ^Tr [dMdM + aU'^ + + ....] , (3.1) 

where M is a x Hermitian matrix. Only the single-trace terms are present in the 
unperturbed theory. We use the standard double-line notations for the propagators and 
vertices. In the t' Hooft limit, g'^N — )• A, all connected correlation functions are of the 
order N"^ with planar diagrams being the leading contributions. Non-planar diagrams are 
suppressed by additional factors of A^. This is easily understood if one recall that each 
vertex contributes a factor of 5"^, each propagator contributes a factor of , and each closed 
loop contributes a factor of A^. Every connected diagram carries an overall factor N^{g'^N)'^ , 
where x is the Euler number for the surface on which the diagram can be drawn. Disconnected 
diagrams that contain multiple connected components carry a factor of A^ that is the product 
of the factors of A'^ for every individual connected piece, and therefore, is larger than A^'^ 
generally. 

Without introducing the replicas, let us first briefly explain what happens if the above 
theory is deformed by a multi-trace term 

5C = ^W[0], W[0] = . (3.2) 

In particular, the double-trace deformation is given by the choice of W[x] = x'^/2. We have 
chosen to normalize the singe-trace operator O slightly unconventionally but compensated for 
it by inserting the additional factor A^^"™ in the multi-trace term so the theory is the same 
to those considered in [9]. The motivation for this choice of normalization is that it reads 
for quadratic deformation, corresponding to a direct coupling between the disorder V 
and the operator O, when interpreted in the context of the replica method. We may choose 
to couple V to some single-trace operator with a different normalization, leading to weaker 
or stronger effect as we will discuss later. 

To help with our counting of the powers of A^ when multi-trace vertices are present, we 
introduce a small trick. Let us split the multi-trace vertex slightly so that it appears as 
separate single-trace ones, as illustrated in Fig. 1, and think of the original vertex as the 
coincide limit of them, i.e. we would take the limit — >• y, z = 1, 2, . . . m, as shown in Fig. 1, 
so they "recombine" into a single multi-trace vertex. This helps because when the multi-trace 
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Figure 1. Momentarily, tiiink of the multi-trace vertex as a group of single-trace ones coincide. 



vertices are involved, strictly speaking one may no longer talk about planar diagrams as the 
vertices themselves can not be drawn on a single plane. Once we split them into multiple 
single-trace ones, the more familiar counting methods can be applied, with only the caution 
that the group of vertices obtained by splitting a multi-trace one do not carry the usual factor 
of g^'^^s and one must be careful with how they are normalized. 

Admittedly, taking the coincide limit of several single-trace vertices to form a multi- 
trace one is not trivial (if one actually does it literally), since various divergences would arise 
in the coincide limit as usual. However, at least within this section, such a splitting and 
recombining process is merely a bookkeeping tool to help us with the counting. One could, 
in principle, rephrase everything below without making the reference to this trick at all, but 
the words become more clumsy as one starts to implement the replica method. Having said 
that. We assume a properly regularization scheme can be applied so that the limit-taking 
process actually makes sense. Since the regularization scheme itself does not involve the 
overall factors of A^, the counting below is not affected. 

Let us consider the correction to the (O) due to a deformation term (3.2) in this theory. 
We treat the deformation perturbatively, so we expand exp{5C} order by order in /, insert 
each term in the path-integral, and evaluate the path-integral with the Lagrangian Eq. (3.1) 
without deformations. To the leading order of /, the correction to ( O ) is given by the integral 

S{0{x)) 

= J d% { 0{x)0{yr r = J d% [{ 0{x)0{yr ) - { 0{x) ) { 0{yy^ )] 

= ^-^^ / lim [( 0{x)0{yi)0{y2) . . . O(y^) ) - ( 0(x) ) ( 0(2/1)0(^2) • • • 0{y^ )] . 

(3.3) 

In the last step, the limit yt y is just a formal notation, but, with those vertices separated, 
we could visualize this expression more easily as Feynman diagrams that connect the vertex 
at X to at least one of those at yi, since diagrams that does not connect x to any of the 2/i's 
are subtracted out. Notice that this is the only requirement for the above expression to be 
non-vanishing and no more vertices are obligated to get connected. Given the constraint, the 
more connected components there are, the more powers of N this correlation function obtains. 
So, the most dominant diagrams that contribute to S {O) must be those that connect only 
X to one of the 2/i's, but leaving all the rest disconnected, as illustrated in Fig. 2. Should all 
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the vertices at y^'s be the ordinary single-trace ones, this diagram would be proportional to 
A^^™ at most. But here, since the group of y^'s together contributes a factor of /A^^"™ in 
stead of A^*", such diagrams are of the order yjv^™~™'+(^^'") = fN'^, exactly the same as the 
usual planar ones. 

We can express the diagrams in Fig. 2, after some proper regularizations, as 

HO) = '^l d%{0{x)0{y)r{Or~' 

where the extra factor of m arises due to the fact that one may choose to connect the vertex 
at X to any one of the vertices at yi leading to equivalent terms. If we translate this object 
by the AdS/CFT dictionary, after proper normalizations, this is the same as 



since 6(3{x)/da{y) is identified with the two-point correlation function {0{x)0{y) Y and /3 = 
( O ). At the same time, the right-hand side of this equation is precisely the change of (3, to the 
leading order of /, due to the change of the boundary condition from a = to a = m//?™"^, 
because 6f3{x) = J d'^yj^^6a{y), and 5a{y) = mf I3{y)"^^^ . The analysis generalizes to 
every higher order corrections similarly, reproducing exactly the same results in the leading 
order of as those given by the altered boundary condition prescription in the holographic 
method. While the analysis here is perturbative in terms of /, the same can be derived 
non-perturbatively by a saddle point calculation [9]. 

We are now at the position to reproduce what has been found in the previous section 
by implementing the replica method in the matrix model, and examine what happens to the 
class of diagrams illustrated in Fig. 2 and why their effects disappear in the limit n — t- 0. 

Let us introduce the replicas so the Lagrangian (3.1) is duplicated by n times, and consider 
order by order the correction to the correlation functions due to the multi-trace deformation 

« r n 

S/:=^-W . (3.4) 
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In the replica method, W[x] = 

Consider S{0), which, to the first order of /, is given by the integral 



6{0{x)) = 




We remind the readers that the correlation functions on the right-hand side are evaluated 
in the theory whose Lagrangian is simply n copies of Eq. (3.1) without the double-trace 
deformation (3.4). Before taking the limit n — )• 0, the disconnected correlation function in 
the last line of the above equation contains several categories of contributions. Feynman 
diagrams that leave all x, yi, and 2/2 disconnected would acquire a factor of n^, diagrams 
similar to those in Fig. 2 that connect only x to one of the yi or 1/2 acquire a factor of n, and 
lastly, diagrams that connect the vertices at all three locations does not acquire any factor of 
n. This is shown expicitly in Appendix B. So, in the limit n — )• 0, only the diagrams in the 
last category survive. 

There is an elegant way to make this discussion clearer and more readily generalizable to 
higher order corrections. Let us consider the following object: 

C = lim - / d^ lim / 0^{x) O.(yi) ^ • • • J] Oiiy^) ) . 

\ i 3 k I I 

As explained in Sec. 2.1, the power of the replica method is that in the limit n — )• 0, it turns 
the above disconnected correlation function into a normalized and fully connected one, i.e. 
by the identity (2.2), we find 

C= lim - / lim ( VOi(x) VO,,(yi) V0fc(y2)--- VOKy,n) 
n-s>o n J yi^y \ ^-^ ^-^ ^-^ ^-^ 

\ i j k I 




hm {0{x)0{yi)0{y2)...0{ym)r • 

y-i-^y 



Here, we arrive at the essential difference between C and Eq. (3.3). The correlation function 
in the last expression above is fully connected before the limit yi — t- y. One may visualize it 
as Feynman diagrams that contain a single connected component even when the multi-trace 
vertex is split. We arrived at this formula quickly by citing the general results (2.2). There, 
it seems essential that all vertices are located at separate positions and the validity of its 
application here might appear questionable. An explicit calculation that makes no reference 
to the splitting and recombining process is presented in Appendix B, which both verifies the 
above formula for the case of a double-trace insertion and serves to clarify what we mean 
exactly by "the fully connected diagrams" here. We remind the readers that in the limit 
yi — >■ y, diagrams shown in Fig. 2 are connected too, but become disconnected once the 
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multi-trace vertex is split. The vitally important difference between those and the correlation 
functions appearing in the last line in (3.5) is that the latter stays fully connected even when 
the multi-trace vertices at y is split into separated single-trace ones! We refer those them as 
the "super connected" diagrams, as illustrated in Fig. 3. 

Should the vertices at yi's be ordinary single-trace ones, the "super connected" diagrams 
would be 0{N'^) because one can draw a planar diagram that connects all vertices. Here, 
however, the counting is somewhat different since the group of vertices at yi contributes 
less powers of N. We already know that diagrams illustrated in Fig. 2 are 0{N'^), which 
immediately implies that "super connected' diagrams being considered here are sub-dominant 
simply because they contain fewer number of connected components. Indeed, in the case of a 
double-trace insertion, a quick counting shows that the leading contributions for d{0) is of 
the order A'^*^. Similar analysis may be carried out to higher order corrections in terms of / 
as well to multi-point correlation functions. 



At this point, it is not surprising at all that if one follows the standard holographic 
prescription for CFTs with double-trace deformations, which takes into account only its 
leading effects at the order of A^^ , one finds all connected correlation functions are not affected 
by disorders once the limit n — )• is taken. In other words, in any CFT that has a well-defined 
large- A'^ expansion such that its correlation functions factorize in the similar way as those in 
the matrix model presented here, the disorder, as it is normalized so far, is just too weak to 
cause any finite corrections to connected correlation functions in the leading order of N. This 
analysis demonstrates that the procedure of implementing the replica method in AdS/CFT 
as we explained in Sec. 2.2 were in fact working properly. 

4 The sub-leading contributions and a new holographic replica method? 

The field theory analysis not only demonstrates that the holographic method gives rise to the 
correct results, but also provides us a mean to look for the sub-leading effects. 

Let us study the order parameter, (O), in the matrix model again. We explained in 
the previous section that, at the order of A^^, the corrections to it due to the double-trace 




Figure 3. The leading corrections to ( O ) when n 0. 
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deformation vanish identically if we take the limit n — )• 0. To the first order of /, the leading 
non-zero effects, described by the class of diagrams shown in Fig. 3, can be written as 

S {0{x)) = I [ d^y hm (0(x)0(yi)0(y2)r ■ 

We had emphasized that splitting the double-trace vertex and then taking the limit that they 
recombine was merely a formal step, which allowed us to easily visualize the difference between 
the diagrams illustrated in Fig. 2 and Fig. 3 and understand why the latter are relatively 
suppressed. The explicit calculation in Appendix B presents an example where the exact 
meaning of such a limit is given. However, if we do take the above expression seriously, we 
may translate such correlation functions back to the AdS side using the standard dictionary. 
The result might be depicted by the Witten diagram shown in Fig. 4 (a). The vertex on the 




(a) (b) (c) 

Figure 4. Corrections to (O) to the first order of /. The vertex y is integrated across the AdS 
boundary. 



boundary at y should be understood as the limit of two vertices yi, y2 — ^ which explains the 
two boundary-to-bulk correlators attached there, and is integrated across the AdS boundary. 
Formally, one may think of it as an insertion of a double trace operator on the CFT side. 
The circular blob in the bulk represents any possible diagrams that connect the two sides. 

Similarly, one could also consider the higher order corrections in terms of /. Again, let 
us take the limit — )• y literally and translate the "super connected" correlation functions 
to the AdS side, we find that one only needs to insert more "double-vertices" on the AdS 
boundary and integrate them all across it, as illustrated in Fig. 5. All such higher order 
corrections are further suppressed by additional factors of N"^ though, so it is not quite useful 
to evaluate them if one does not at the same time include all other competing contributions 
such as the loop diagrams in the AdS bulk. 

In passing, we mention an observation that follows immediately from the above analysis. 
When the bulk Lagrangian is described by Eq. (2.9) and the disorder is coupled to the 
operator O, to the leading order of /, the correction to the condensate ( O ) only shows up if 
the system is in the superconducting phase, i.e. ( O ) 7^ 0, a result fully expected of course. 
This is because the Witten diagram in Fig. 4 (a) can not be completed without breaking the 
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(a) (b) 
Figure 5. Higher order corrections to (O). Vertices yi,2,3 are integrated across the AdS boundary. 

U{1) gauge symmetry. Only if ip has a non-zero background profile, such diagram can be 
realized, as shown in Fig. 4 (b), if one of the legs is connected to the background. Fig. 4 (c) 
gives a more explicit realization. All those diagrams are accompanied by their conjugate of 
course. 

We could also compute the corrections to the two-point correlation functions perturba- 
tively in / in the same way. For example, the leading corrections to the gauge field boundary- 
to-boundary propagator is given by the diagrams shown in Fig. 6. 




(a) (b) (c) 

Figure 6. The corrections to the averaged boundary-to-boundary correlator for the gauge field to the 
first order of /. The vertex z is integrated across the AdS boundary. 

The methods proposed here deserve some discussions. Treating every insertion of a 
double-trace operator on the CFT side by simply adding a "double vertex" on the AdS 
boundary for the dual field seems superficial. If not in the context of the replica method, 
such a map is mistaken. After all, whenever the number of replicas is greater than zero, 
the leading contributions are given by the diagrams show in Fig. 2, not at all included in 
the formalism suggested here. However, if one is only interested in the limit n — t- 0, the 
missing contributions all disappear, and it becomes possible to described all the surviving 
contributions by the simple Witten diagrams mentioned. 

Some of the diagrams shown in Fig. 4 are similar to the I-loop diagrams considered in 
[11], but milder since only one of the vertices of the "loop" needs to be integrated throughout 
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the AdS bulk. They will not suffer from the extra "volume" divergences discussed in [11] for 
this reason. It is only the d-dimensional Minkovski part of the loop integrals that may need 
regularizations and the standard QFT methods in flat spacetime should suffice. 

Apart from possible regularizations needed, the true weakness of the procedure outlined 
here is that all higher order corrections are suppressed by more powers of and it would 
be pointless to evaluate them if one does not include all the other competing effects at the 
same time. Let us mention, however, a potentially much more interesting scenario. It is 
conceivable to make the effects of the disorder stronger so they show up at the leading order 
of the large-A expansion, i.e. 0{N'^). To do so, we just need to make the disorder stronger 
as A — 7- oo, leading to a double-trace deformation with a different normalization: 

5C=^0'^. (4.1) 

Formally, we can achieve this goal in the matrix model if we just let the coupling between V 
and O depend on the gauge coupling g. More specifically, we may put the same overall factor 
of g^"^ in front of the coupling of V and O so it reads 

5Lv = ^V{x)0{x). 

Integrating out V leads to a deformation term 5C = tJ-^O"^ and in the t' Hooft limit it 

becomes 5C = ^^jrO"^- After redefining A~^/ — )• /, we find (4.1). Equivalent, we may 
change the Gaussian distribution profile for V so that the width becomes /A^. In any case, 
it amounts to choose a different normalization for the single-trace operator coupled to the 
disorder. By giving / an extra factor A^, the leading diagrams shown in Fig. 3 are made 
0(A2). 

Normally, a double-trace term normalized as in Eq. (4.1) is disastrous, because it destroys 
the large-A expansion. As we bring up the sub-dominant contributions to the leading order, 
we necessarily make all the leading ones, shown in Fig. 2, too large simultaneously. As A — t- 
oo, they totally over-dominate any other physics. If one consider the higher order corrections 
in terms of /, situation only becomes worse and worse, essentially making the large-A limit 
singular. Consequently, one does not expect to find a holographic dual describable by a 
semi-classical field theory. 

In the current context, however, we are rescued by the magic power of the replica method. 
Because we are only interested in the limit n — t- and those bad-behaved diagrams are known 
to vanish in this limit, it is hopeful that the said disease can be cured. Indeed, the diagrams 
as shown in Fig. 2 are at best of the order nN^ as one can easily verify in the example given 
in Appendix B. Therefore, it suffices if we somehow restrict ourselves in the corner of the 
parameter space, such that, as A^ — )• oo, n = nN"^ is kept finite. That way, the diagrams of 
the order nA^ become 0{N'^) as well. It is not hard to check the pattern is maintained when 
higher order corrections in terms of / are considered. Setting n = leaves with us only the 
diagrams similar to those illustrated in Fig. 3, but leaving h ~ 0(1) allows both diagrams in 
Fig. 2 and Fig. 3 to contribute at the same order. 
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Interestingly, with the normahzation in Eq. (4.1), it also happens automatically that the 
higher order contributions in terms of the deformation strength / are no longer relatively 
suppressed. All of them become 0{N'^) and contribute at the leading order in the large- 
N expansion simultaneously, allowing themselves to stand out and be separated from other 
competing effects, such as those of the non-planar diagrams, or correspondingly loops in the 
AdS bulk. In this theory, both / and n can be considered as an additional 't Hooft coupling. 

This observation also implies that it becomes both necessary and meaningful to sum all 
the contributions of the double-trace deformation to all orders of /, which is out of our reach 
at the moment unfortunately. It is plausible that a modified AdS/CFT correspondence, yet 
unknown, may exist that capture all the contributions non-perturbatively. Somehow the dual 
theory must incorporate n in a non-trivial way because only when it is kept 0(1), a consistent 
large-A^ expansion can be defined. 

5 Discussion and Conclusions 

While we used the matrix theory as a toy model to probe the effect of multi-trace deformations 
perturbatively on the field theory side, the discussion in the previous sections solely relies on 
the essential property that the correlation functions factorize in the large- expansion. We 
believe that the found sub-leading contributions remain the same for general theories with 
holographic duals. 

For realistic physical questions, however, one must restrict himself to the time-independent 
random disorders, which is equivalent to giving the disorder strength / a factor of d{oj) where 
uj is the frequency of the disorder. For numerical calculations, one must then smear out 
this delta-function somewhat, replace it by a sharp Gaussian-like profile, and take the limit 
properly in the end. 

In our discussions above, we had not specified the spacetime signature. What we found 
obviously applies in the Euclidean space where the correlation functions are defined as written. 
In the real-time formalism, however, various type of correlation functions can be defined, and 
physically the retarded ones are the most interesting. Our discussions above should apply to 
all correlation functions, since the counting of the powers of A^ is a separate issue. On the field 
theory side, different correlation functions correspond to different choice of the time-contour 
in the path-integral, and, on the AdS side, the choice amounts to specifying different boundary 
conditions for the mode solutions at the horizon [28-31], and the discussions presented above 
remain unchanged. 

It is a bit disappointing that the method proposed here for evaluating the non-vanishing 
sub- leading effect is only perturbative in terms of the disorder strength /. In fact, it is simply 
pointless to evaluate the higher order corrections on their own, since being more and more 
suppressed in the large- A^ expansion, they are contaminated by more and more competing 
contributions not discussed in this paper. 

A much more interesting scenario is found if we decide to increase the disorder strength 
such that, in the replica method, it leads to a double-trace deformation normalized as in 
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Eq. (4.1). In an ordinary CFT, such a normalization is pathological, since connected diagrams 
in the deformed theory may carry arbitrarily higher powers of A^, totally invalidating the 
large-A^ expansion. However, as we put the formalism in the context of the replica method, 
the situation becomes very different, because there emerges a new parameter that is small, 
namely, the number of the replicas. We found that the problematic diagrams that may carry 
higher powers of N must always be proportional to {nN'^)'^N'^. So, focusing on the limit 
n — 7- allows us to consider this theory more seriously. In fact, as long as we can keep 
h = nN'^ finite as — t- oo, all the leading connected diagrams are 0{N'^) and a consistent 
large- expansion can be defined, in which n may be thought of as a new 't Hooft coupling. 
In the meanwhile, higher order contributions in terms of / are also no longer suppressed, and 
they all contribute at the order oiN"^. Therefore, at every order of /, the leading contributions 
due to the multi-trace operator can be cleanly separated from other corrections such as those 
given by the non-planar diagrams. It would be extremely interesting, if one could find a new 
holographic model for such a theory that captures all the effects non-perturbatively. This is, 
of course, a novel regime since the duality may only work properly if we can keep the number 
of the fields, n, to be 0{N~'^) as — )• oo, something that remains unclear to us how to 
achieve but is surely worth further investigations. 
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A The replica method for free theories 

The simplest example that the replica method can be explicitly implemented is a free theory 
where the operator coupled to the disorder is the fundamental field itself. It is unnecessary 
to elaborate such a simple example if not because we find that the holographic calculation 
surprisingly resembles the free theory in various ways. 

Let us consider the action 5o = —\KX'^^ where we are being schematic by suppressing 
space-time integrals and denote the kinetic operator simply by K. Let it be coupled to 
a random potential V so we have S = Sq + VX and let V has a distribution given by 
P\V] = exp{ — ^y^}. By the replica method, we should consider the following action: 
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To deal with the mixing term, we define Xi = CLijXj, where aij is the matrix defined in 
Sec. 2.1. In terms of Xi the action is diagonafized and we have 



S = -lEKXi + ^Xl 



i 



So we easily find 



(X) = lim-( VX, ) = Iim/X„\ =0. 




{XX) =lim- [(l„X„) + (n-l)(l„_iX„_i)] =- + -A, 



{X){X)— (^XnXn'j — (^Xn-lXn-l'j ~ J{2 ^ 



but the connected correlation function, given by the difference of the two expressions above, 
is independent of / and 



This is not at all surprising. Since the action is quadratic, the connected two-point function 
for X is independent to the linear term XV, so for any chosen V, the correlator remains to 
be 1/K, and no wonder the averaged value is the same. The essential point that makes such 
simple calculations possible is that the system can be easily diagonalized by the rotation aij. 
We see that the same holds true for correlation functions in the context of the holographic 
replica method. 

B An explicit calculation that verifies Eq. 3.5 

While Eq. (3.5) is general, it is not fully self-evident. We present an explicit calculation in 
the case of a single double-trace insertion, which not only verifies the Eq. (3.5) but also helps 
to clarify the actual meaning of the final expression given there. 

Let us consider, for example, the leading order correction to ( O ) given by the integral: 



Here, all the correlation functions on the right-hand side are evaluated in the limit / = 0, 
and so all the replicas decouple. Clearly, we have 




ma)0,{b)r = Z^ma)){0,{b)) 
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for Vi / j. Here Z = J DX exp{S[x]} is the partition function for the single-copied theory, 
and VEVs on the right-hand side of the above equation are normahzed. So to the leading 
order of /, we find 

<^ ( Oi{x) )° = ^ lim ^ / d% [{ 0,{x)0,{y)Oi{y) ) + (n - 1) ( Oi(x) ) ( 02{y)02{y) ) 
2 n->0 n J 

+ 2(n-l)(Oi(x)Oi(y))(02(y)) 

+(n - l)(n - 2) ( Oi(a;) ) ( 02(y) ) ( Os{y) )] 



d^2/(0(x)0(y)0(2/)r, 
where 

{ABCy = {ABC)-{AB){C)-{A){BC)-{AC){B) + 2{A){B){C) 
is precisely the totally connected correlation functions if A, B and C are all separated. 
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